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Abstract We review hadron formation from a deconfined quark gluon plasma (QGP) via coa- 
lescence or recombination of quarks and gluons. We discuss the abundant experimental evidence 
for coalescence from the Relativistic Heavy Ion Collider (RHIC) and compare the various coa- 
lescence models advocated in the literature. We comment on the underlying assumptions and 
remaining challenges as well as the merits of the models. We conclude with a discussion of some 
recent developments in the field. 

CONTENTS 



Annu. Rev. Nucl. Part. Sci. 2008 72 

Introduction 

Hadronization 

Early Approaches to Recombination 

Challenges at RHIC 



Formulations of Hadronization by Recombination 

Basic Theory 

Different Implementations of Recombination . . . . 
Competing Mechanisms of Hadron Production . . . 

Elliptic Flow 

Comparison of Approximations and Assumptions . 



Challenges and Outlook .... 

Energy and Entropy 

Space- Momentum Correlations 
Outlook 



Data from Elementary and Heavy-ion collisions 24 

Hadron Spectra and Baryon to Meson Ratios 

Elliptic Flow and Quark Number Scaling 

Heavy Quarks 

Particle Correlations and Fluctuations 

Beam Energy Dependence 



Conclusions 37 



Acknowledgments 



2 



Hadronization by Coalescence 3 
1 Introduction 

Collisions between heavy nuclei are used to probe the properties of nuclear mat- 
ter at high temperature and density. Lattice QCD calculations indicate that if 
nuclear matter is heated above a critical temperature Tc ~ 185 MeV, quark and 
gluon degrees of freedom will be liberated and a deconfined quark-gluon plasma 
(QGP) forms ([I112])- Unambiguous signatures of quark gluon plasma formation 
in heavy-ion collisions have been sought for decades. Recently, experiments at 
the Relativistic Heavy Ion Collider (RHIC) have presented evidence that such a 
new state of matter has finally been found in collisions of Au atoms at a center 
of mass energy of ^/s = 200 GeV per nucleon-nucleon pair (OH]). 

The hot QGP phase formed in nuclear collisions at RHIC with a core temper- 
ature in excess of 300 MeV only lasts for an extremely short time. It quickly 
expands due to the large pressure and cools on the way. Eventually, the quark 
and gluon constituents need to combine into color-neutral objects and hadrons 
have to be formed when the temperature reaches Tc. The process of hadroniza- 
tion from a QGP may be quite different from hadronization in other cases, such 
as hadronization of hard scattered parton in elementary collisions where no ther- 
malization is reached and no bulk of partons is formed. In this review we discuss 
a model of QGP hadronization by coalescence or recombination of quarks and 
gluons. The models discussed here have had success in describing many salient 
features of hadron production in heavy- ion collisions. 

The emergence of recombination models was largely motivated by several unex- 
pected observations ([Sj) which were discussed as "the baryon puzzle" for a while. 
This was referring to measurements of baryon production in the intermediate 
transverse momentum region (1.5 < < 5 GeV/c) ([Sill]). Both the yield and 
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the elliptic flow of baryons exhibited strange features. In nucleon-nucleon colli- 
sions at = 3 GeV/c, only one baryon is produced for every three mesons (1:3), 
reflecting the larger mass and the requirement of a non-zero baryon number to 
form the baryon. In Au-|-Au collisions at RHIC however, baryons and mesons are 
created in nearly equal proportion (1:1) despite those differences. In the same 
PT-region, the elliptic anisotropy (f 2) of baryons is also 50% larger than that for 
mesons. Therefore, baryon production is particularly enhanced in the direction 
of the impact vector between the colliding nuclei (in-plane) (HKT]). 

The large baryon V2 eliminates several possible alternative solutions put for- 
ward for the baryon puzzle. The most common explanations for the baryon 
anomaly at RHIC were 

• coalescence or recombination — Multi-quark or gluon processes during 
hadron formation piQlfTOlfTTtn^ . 

• baryon junctions — Gluon conflgurations that carry baryon number (jl3|) . 

• flow — Collective motion that populates the higher pj^-regions of phase 
space for the more massive baryons, as described by hydrodynamics ()141 

Only coalescence models have survived the tests imposed by an impressive amount 
of data taken after the original discovery of the baryon enhancement. They are 
particularly attractive because they seem to provide a natural explanation for the 
valence quark-number scaling that has been observed in V2 measurements. They 
also relate hadronic observables to a pre-hadronic stage of interacting quarks 
and gluons. As such, they touch on questions central to the heavy-ion physics 
program: deconfinement and chiral symmetry restoration. 

This review is organized as follows. In the remainder of this section we discuss 
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the general context of hadronization and a brief history of recombination models. 

We also review the experimental evidence from RHIC. In Sec. [2] we review the 

basic theory and compare the different implementations of recombination models. 

In Sec. [3] we present a comprehensive overview of the available data which can 

be addressed by coalescence. We conclude with a discussion of open questions, 

recent developments and future directions of research in Sec. [H 

1.1 Hadronization 

Hadronization has always been a challenging aspect of quantum chromo dynamics 
(QCD), the fundamental theory of the strong force. QCD bound-states are non- 
perturbative in nature and a first-principle description of their formation has yet 
to be obtained. In this subsection we briefly discuss two approaches to deal with 
hadronization which are routinely used in nuclear and particle physics; both of 
them have connections to the recombination model discussed in this review. 

Light cone wave functions are used to describe the structure of hadrons relevant 
for exclusive processes <\17\i . Exclusive here means that they deal with a full set 
of partons with the quantum numbers of the hadron. Exclusive processes at high 
momentum transfer are naturally dominated by the few lowest Fock states. For- 
mally, light cone wave functions are matrix elements of the set of parton operators 
between the vacuum and the hadron state in the infinite momentum frame, e.g. 
(j)p ~ {0\uud\p), schematically, for a proton p. They describe the decomposition 
of the hadron in longitudinal momentum space in terms of partons with momen- 
tum fractions Xj. From theory these wave functions are only constrained by very 
general arguments like Lorentz-covariance and approximate conformal symmetry. 
Direct measurements are difficult, but estimates have become available in recent 
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years, in particular for the lowest Fock state of the pion (|18 yi9p . 

A complementary technique has been developed for inclusive hadron produc- 
tion, initial state — > h + X, at large momentum transfers in which a single colored 
parton a has to hadronize into the hadron h. For this purpose fragmentation or 
"parton decay" functions Da^h{z) have been defined ([20]) . They give the proba- 
bility to find the hadron h in parton a with a momentum fraction z, < 2 < 1. 
The cross section for inclusive hadron production in + e~, lepton-hadron or 
hadron-hadron collisions can then be written as a convolution 

(TH = CFa® Da^h (1) 

of the production cross section CTa for parton a with the fragmentation function 
Da^h{z) ()2ip22p . Fragmentation functions are not calculable in a reliable way 
from first principles in QCD. However, they are observables and can be measured 
experimentally. Parameterizations using data mostly from e"*" -|- e~ collisions are 
available from several groups (|23p . Physically, the fragmentation of a single 
parton happens through the the creation of qq pairs (through string breaking or 
gluon radiation and splitting) which subsequently arrange into color singlets, and 
eventually form hadrons. 

Both examples above apply to processes with a large momentum transfer, i.e. 
with a perturbative scale ^ S> Aqcd- They are based on the concept of QCD 
factorization which separates the long and short distance dynamicso Such a 
perturbative scale is absent for the hadronizing bulk of partons in a heavy ion 
collision and neither technique, fragmentation nor exclusive wave functions, can 



^Note that we have neglected the scale dependence in the notation for wave functions and 
fragmentation functions for simplicity. A discussion of the scale dependence can be found in the 
original references given. 
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be readily applied in this situation. 

To see the challenge more clearly, let us compare the different initial condi- 
tions for the hadronization process. Fragmentation applies to a single parton in 
the vacuum, whereas exclusive wave functions are applied to a full set of valence 
quarks in the vacuum. On the other hand, the initial state just before hadroniza- 
tion in nuclear collisions is a thermal ensemble of partons just above Tc. The 
exact degree of thermalization is not clear a priori, but we will see below, that 
complete thermalization might not be necessary. 

Rather, the crucial point seems to be that partons have a certain abundance 
in phase space such that there is no need for the creation of additional partons 
through splitting or string breaking. The most naive expectation for such a sce- 
nario is a simple recombination of the deconfined partons into bound states. In- 
deed, there is experimental evidence that this is the correct picture for hadroniza- 
tion even long before a thermal occupation of parton phase space is reached. 

1.2 Early Approaches to Recombination 

Recombination models have first been suggested shortly after the invention of 
QCD in the 1970s. They successfully described hadron production in the very 
forward region of hadronic collisions (24). The observed relative abundances of 
hadrons clearly deviate from expectations from fragmentation in this region. This 
is known as the leading particle effect ()25p . E.g. a clear asymmetry between D~ 
and mesons was found in fixed target experiments with vr" beams on nuclei 
by the FNAL E791 collaboration (26). The measured /D^ asymmetry goes to 
1 in the very forward direction, while fragmentation predicts that this asymmetry 
is very close to 0. This result can be explained by recombination of the c from a cc 
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pair produced in the collision with a d valence quark from the beam vr" remnants. 

This mechanism is enhanced compared to the c+d recombination which involves 

only a sea quark from the 7r~ (|27p . There is no thermalized parton phase in this 

example, which strongly backs our argument at the end of the last subsection. 

We are led to the important conclusion that the presence of any reservoir of 

partons leads to significant changes in hadronization. Vacuum fragmentation is 

no longer a valid picture in this situation. The reservoir of partons in the case 

of the leading particle effect is the soft debris from the broken beam hadron. In 

heavy ion collisions it is the distribution of thermal partons. First applications 

of the coalescence picture to nuclear collisions appeared in the early 1980s (f28|) . 

This eventually led to the development of the ALCOR coalescence model in the 

1990s (j29y 30 t[5T|) . ALCOR focuses on hadron multiplicities and was successfully 

applied to hadron production at RHIC and the lower energies at the CERN SPS. 

1.3 Challenges at RHIC 

Results from the first years of RHIC triggered a revival for recombination models 
applied to heavy ion collisions in an unexpected region. Three measurements in 
particular, taken in the intermediate px range (1.5 GeV/c < pt < ^ GeV/c), 
have defied all other explanations. This region is outside of what was thought 
to be the "bulk" of hadron production (pT < 1.5 GeV/c) whose features should 
be described by thermalization and hydrodynamic collective motion (ALCOR 
describes bulk hadronization). Rather, the intermediate pT region was expected 
to be dominated by fragmentation of QCD jets, after it was confirmed that this 
was the case for pion production in p + p collisions at RHIC (32). However, the 
results from RHIC clearly pointed towards a strong deviation from the fragmen- 
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tation process at intermediate px in central Au+Au collisions. The three key 

observables were 

• the enhanced baryon-to-meson ratios ([5l [33|) . 

• the nuclear modification factors Raa and Rcp — i-e. the ratio of yields 
in central Au+Au collisions compared to peripheral Au+Au {Rcp) or 
p + p (Raa) collisions scaled by the number of binary nucleon-nucleon col- 
lisions (5,6). 

• the anisotropy of particle production in azimuthal angle relative to the 
reaction plane — i.e. the elliptic flow parameter V2 <!^\S5 \ \3 ^ \7]l . 

Fig. [1] shows the measured anti-proton/pion {5) and ^/K^ ([33|) ratios as a 
function of pT for various centralities and collision systems. At intermediate pT, 
a striking difference is observed between the baryon-to-meson ratios in central 
Au+Au collisions and those in e'^ + (36D or p + p collisions (|37p. The mea- 
surements in Fig. [T] indicate that the process by which partons are mapped onto 
hadrons are different in Au+Au collisions and in p + p collisions. Changes solely 
to the parton distributions prior to hadronization are not likely to lead to such 
drastic changes in the relative abundances. 

Fig. [2] shows the nuclear modification factor Rcp measured at RHIC for various 
identified hadrons. If the centrality dependence of particle yields scales with the 
number of binary nucleon-nucleon collisions, Rcp will equal one. A suppression 
at high pt is taken as a signature for the quenching of jets in the bulk matter 
formed in central collisions. However, baryons (A + A, S + H, and 17 + fi) (I33y38p 
systematically show less suppression than mesons (kaons or (p) ()39tl40p ). The 
same behavior was found for protons and pions (j^T]) . This key result shows that 
the mass of a hadron is less important for its behavior at intermediate pT than 
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the fact whether it has two or three valence quarks. This ruled out explanations 

blaming collective motion (flow) for the baryon enhancement, and it is a strong 

indication that parton degrees of freedom are important. Last doubts were erased 

by a direct comparison of protons and (p mesons which have the same mass but a 

different valence quark content: (p mesons behave like other, lighter mesons, not 

like protons (Ii0l[i2l) . 

In non-central nucleus- nucleus collisions, the overlap region of the nuclei is el- 
liptic in shape. Secondary interactions can convert this initial coordinate-space 
anisotropy into an azimuthal anisotropy of the final momentum-space distribu- 
tion. That anisotropy is commonly expressed in terms of the coefficients from 
a Fourier expansion of the azimuthal dependence of the invariant yield (j43p . see 
Sec. [21 The second component (the "elliptic flow" parameter V2) is large due to 
the elliptic shape of the overlap region. Fig. [3] shows the measured values for V2 
as a function px for pions, kaons, protons and Lambda hyperons dHUSSD- In the 
bulk region (px < 1-5 GeV/c), V2 is increasing with pT (jSj). In this region the V2 
values for different hadrons are ordered by their mass with more massive particles 
having smaller V2 values (j6ll35p45p . This mass ordering is qualitatively understood 
in hydrodynamic models of the expansion of the bulk of the fireball |14j). Some 
hydrodynamic calculations are also shown in the figure. For px > 1.5 GeV/c, the 
data clearly deviates from hydrodynamic calculations. The measured V2 seems to 
saturate, as predicted by parton cascades ((Mj), and the particle- type dependence 
reverses: V2 values for the more massive baryons are larger than those for mesons. 

V2 can also be generated if jets are quenched in the quark gluon plasma (|T7l) . 
However, such calculations grossly underestimate the measured values of V2, in 
particular when they simultaneously have to explain values of Rcp close to one. 
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The data clearly shows that although protons and hyperons have Rcp values near 
unity, their maximum V2 values exceed those of pions and kaons by approximately 
50%. Taken together, the particle-type dependence of V2 and Rcp provide very 
stringent tests of various models for particle production and have ruled out pure 
jet fragmentation or simple hydrodynamics as models for hadron production at 
intermediate pT- 

2 Formulations of Hadronization by Recombination 
2.1 Basic Theory 

Coalescence or recombination of particles is a very general process that occurs in 
a wide array of systems from the femtometer scale to astrophysics. In all these 
fields a first approach is to discard the details of the dynamical process in favor 
of exploiting an adiabatic approximation in which a projection of the initial state 
onto the final clusterized state is considered. In the specific case of recombination 
of partons, most work found in the literature uses an instantaneous projection 
of parton states onto hadron states. The expected number of hadrons h from a 
partonic system characterized by a density matrix p is given by 



Instantaneous here means that the states are defined on a hypersurface which 
is typically either taken to be at constant time, t = const., or on the light-cone 
t = ±z. In this case information about the hadron bound state is schematically 
encoded in a wave function or Wigner function. As we will see this approach 
leads to very simple math, but it has the conceptual disadvantage that only 
three components of the four momentum are conserved in such a2— >lor3— ^1 
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coalescence process. A more dynamical approach based on resonance scattering 

can be realized, which avoids this problem (|48p. The information about the 

hadron bound state is then encoded in a cross section. In this section, we will 

focus on the instantaneous projection formalism which has had great success 

explaining RHIC data. We will come back to the dynamic formulation in Sec. [H 

All available models of instantaneous coalescence can be traced back to the 

following basic formula which can be derived from Eq. ([2]). The number of mesons 

with a certain momentum P is (1491) 



dNM V- f d^R d?qd?r ( r P r P \ .... 

15F = 5 / W'*-- - 2' 2 - + r T + <l) *"'(^. ■!>• (3) 

Here M denotes the meson and a, h are its coalescing valence partons. Wah 

and $jv/ are the Wigner functions of the partons and the meson respectively, P 

and R are the momentum and spatial coordinate of the meson, and q and r 

are related to the relative momentum and position of the quarks. The sum runs 

over all possible combinations of quantum numbers of the quarks in the hadron, 

essentially leading to a degeneracy factor Cm- 

Note that coalescence, just as its counterpart in exclusive processes, is based 
on the assumption of valence quark dominance, i.e. the lowest Fock states are 
the most important ones. The corresponding formula for baryons containing 3 
valence quarks can easily be written down as well. It is also straightforward to 
generalize Eq. ([3]) to include more partons which would be gluons or pairs of sea 
quarks, accounting for the next terms in a Fock expansion ()5ip . 

For a meson consisting of two quarks its Wigner function is formally defined 

as 

$M(r, q) = I d^.e— Va/ (r + |) (r - (4) 
where the 2-quark meson wave function in position space (/9m can be represented 
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as 

(n; r2|M; P) = e-'^^'^+^'^^^M (ri - ra) (5) 

The Wigner function of the partons can be defined in a similar way from the 
density matrix p (j49p . 

To evaluate Eq. ([3]) expressions for the hadron wave functions and for the 
distribution of partons have to be used as input. We discuss the different imple- 
mentations in the next subsection. Let us emphasize two common features of all 
implementations. For one, the Wigner function for the multi-parton distribution 
is usually approximated by its classical counterpart, the phase space distribution 
of the partons on the hypersurface of hadronization. Secondly, Eq. ([3]) is made 
explicitly Lorentz-covariant to account for the relativistic kinematics. 

2.2 Different Implementations of Recombination 

Different manifestations of Eq. ([3]) have been used in the literature (]52p. Clos- 
est to the master formula is the implementation by Greco, Ko and Levai [GKL] 
P HI53p . In this approach the full overlap integral in Eq. ^ over both relative 
position and momentum of the partons is calculated. On the other hand, sev- 
eral groups, (e.g. Fries, Miiller, Nonaka and Bass [FMNB] ([inilMlEO]) ; Hwa and 
Yang [HY] ([5 ^155]) and Rapp and Shuryak [RS] (56) simplify the situation by 
integrating out the information about position space. This leads to a formula- 
tion solely in momentum space in which the information about the hadron is 
further compressed into a squared (momentum space) wave function (also called 
a recombination function by some authors). 

The implementation by Greco, Ko and Levai was originally motivated by a rel- 
ativistic extension of the formalism for the coalescence of nucleons into deuterons 
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and other light clusters in relativistic heavy-ion collisions. Coalescence has been 

successfully applied to nucleons for more than two decades ([57 1158]) . GKL use a 

manifestly covariant version of Eq. ([3]) for the number of mesons coalescing 

Nm = Cm J Yl(P' d'^Pi ^(pI ~ mlWab{ra,Pb; rb,Pb)^M{r; q) ■ (6) 

i=a,b 

The relative phase space coordinates r = ri, — Va and q = Pb — Pa are the four- 
vector versions of the vectors r and q in Eq. ([3]). da is a volume element of a 
space-like hypersurface. The hypersurface of coalescing partons is usually fixed 
by GKL through the condition of equal longitudinal proper time r = Vf^ — z^- 

In the GKL formalism the full phase space overlap of the coalescing particles 
is calculated. For mesons this leads to a 6-dimensional phase space integral 
which is computed using Monte-Carlo techniques (|53p . This has the advantage 
to avoid some of the more restrictive approximations employed by other groups. 
In addition, the numerical implementation of the 6D-phase space integral can 
be applied directly to a quark phase which has been extracted from a realistic 
dynamic modeling of the phase space evolution in the collision. Soon after the 
first implementation of GKL, similar techniques were used for hadronization in 
the partonic cascade approach by Molnar (|59p . 

The hadron Wigner function for light quarks used by GKL is a simple product 
of spheres in position and momentum space 

Qtt r 1 r 1 

^A/(r; q) = ye - r2] X G [A^ - + (mi - ^2)^] . (7) 

The radii and Ap in the Wigner formalism obey the relation Ap = A~^, 
motivated by the uncertainty principle. The parameter Ap is taken to be different 
for baryons and mesons and is of the order of the Fermi momentum. Gaussian 
Wigner functions are used if heavy quarks are involved (|60l) (see also Fig. [9]) . 
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The Wigner functions used in the GKL approach appear to be more arbitrary 

then those based on hght cone wave functions (see FMNB below) However, it 
turns out that many aspects of coalescence do not depend critically on the wave 
function for systems close to thermalization. On the other hand, using the full 
information about the phase space distribution of partons permits a direct connec- 
tion to many quantitative properties of the bulk of the fireball, like the multiplicity 
of partons, and the energy and entropy densities just before hadronization. The 
parameters found by GKL in order to reproduce the behavior of hadron spectra 
at intermediate pT also provide a bulk of the partonic fireball which is consistent 
with what can be inferred from hydrodynamics and experimental data. E.g., the 
radial flow, parameterized as /? = Por/R exhibits a slope parameter /3o = 0.5 con- 
sistent with hydro dynamical calculations at the end of the quark-gluon plasma 
phase (|15|) . Moreover the energy density at hadronization is e = 0.8 GeV/fm^, 
which is very close to what is expected from lattice QCD calculations (PQE]). In 
addition, the entropy is found to be dS/dy ~ 4800 in agreement with the value 
inferred from experimental data by Pratt and Pal (j6ip . 

Compared to the full phase space implementation which is rather complex, 
simplified momentum space models focus on a direct exposition of some key fea- 
tures which can then be treated analytically. We will discuss the implementation 
by Pries, Miiller, Nonaka and Bass in detail here, but the approaches taken by 
Hwa and Yang, and Rapp and Shuryak are very similar. 

The first assumption made by FMNB is that variations in the quark distribu- 
tion across the size of a pre-hadronic state (which may be smaller than a free 
hadron) are small. The integration over the relative position of the quarks can 
then be carried out. For further simplification one focuses on the case that the 
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momentum |P| of the hadron is much larger than the mass M. This ahows 
one to treat the hadron as being on the Ught cone with a large + -momentum, 
P+ 3> P~ (the z-axis in the lab frame is here taken to point into the direction of 
the hadron; this is called the hadron light cone frame in (J49')). The momenta of 
the partons inside the hadron can be parameterized by light cone fractions Xi of 
(0 < Xj < 1) and transverse momenta kj orthogonal to the hadron momen- 
tum P. The momenta kj are usually integrated as well in a trivial way, leaving 
a single longitudinal momentum integration. 

In the absence of any perturbative scale the light cone wave functions are not 
known from first principles. But, again, the coalescence from partons thermally 
distributed in phase space is not very sensitive to the shape of the wave functions. 
For the lowest Fock state of a meson the squared wave function or recombination 
function is usually parameterized as (j49|) 



Here the Ui are powers which determine the shape, and the constant B is fixed 
to normalize the integral over $ m to unity. The yield of mesons with momentum 
P can then be expressed as 



where da is the hypersurface of hadronization. In many cases the emission integral 
over the hypersurface is not calculated explicitly, but replaced by a normalization 
factor proportional to the volume of the hadronization hypersurface. 

Several choices for the powers can be found in the literature. Asymptotic 
light cone distribution amplitudes suggest = 2 for light valence quarks. For 
heavy-light mesons the relative size of the powers has to be adjusted such that the 




(8) 
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average velocity of the quarks is about the same. E.g. for a charm and hght quark 

system, like the D meson, values Oc = 5, Uu^d = 1 are used by Rapp and Shuryak 

(|56p . It is sometimes useful to look at the extreme case a, ^ oo with the ratio of 

the Oi fixed. For two light quarks this implies ^Mi^i, X2) = 6{xi — l/2)6{x2 — l/2). 

This is the limit of a very narrow wave function in momentum space and the 

remaining integral is trivial. 

Analytic implementations of recombination applied to intermediate pT in heavy 

ion collisions usually assume a thermal distribution of partons at hadronization. 

For such a system it seems to be sufficient to neglect correlations between partons 

and to use a factorization into single-particle phase space distributions 

Wab{ra,Pa]rb,Pb) = fa{ra, Pa) fbin, Pb) ■ (10) 

With thermal one-particle distributions / this gives good results for single inclu- 
sive spectra, hadron ratios, etc. at RHIC. Observables dealing with correlations 
of hadrons are more sensitive to correlations among partons. Results from RHIC 
seem to suggest that jet-like correlations between bulk partons exist down to in- 
termediate pt and have to be taken into account (62). This will be addressed in 
more detail below. Coalescence applied to non-thermal systems, e.g. to parton 
showers ()63t[M|) . or hadronic collisions ([56|) . require more sophisticated models 
for multi-parton distributions. 

As we have mentioned above, recombination of partons in a thermal ensemble 
has the interesting property that the process is largely independent of the shape 
of the hadron wave function. This can be most easily seen using the FMNB 
formalism with the partons coming from the tail of a Boltzmann distribution 
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e For a meson with large momentum P the integral in ([9|) is 



1 



1 



dXadXh^M{Xa,Xb)e 



XbP/T 




dXadXb^M{Xa,Xb) , (H) 



Jo 



independent of the shape of ^m- Moving to lower hadron Pt or using the full 
phase-space overlap, as in GKL, make this argument less rigorous. But even 
in those cases the results are only weakly dependent of the shape of the wave 
function unless very extreme choices are made. From this small exercise we 
can read off another important fact, which is tantamount to solving the baryon 
puzzle at RHIC: both mesons and baryons would lead to the same Boltzmann 
distribution ~ e~^/"^ (for sufficiently large momentum P), which is very different 
from the suppression for baryons expected from fragmentation. 

2.3 Competing Mechanisms of Hadron Production 

In order to compute realistic hadron spectra that can be compared to data mea- 
sured in heavy-ion collisions, other important mechanisms of hadron production 
have to be considered as well. QCD factorization theorems state that leading- 
twist hard parton scattering with fragmentation is the dominant mechanism of 
hadron production at asymptotically high momentum transfer (22). This can 
also be seen from the simple analytic formulas discussed in the previous subsec- 
tion. Let us again consider the tail of a thermal parton distribution /th ~ Ae~^/^ 
and compare it to a power-law distribution /jet ~ Bp~°' for large p. Power-law 
distributions are typical for partons coming from single hard scatterings. Both 
recombination and fragmentation preserve the basic shapes of the underlying 
parton distribution. 

As we have already argued above, recombination of n thermal partons leads to 
a thermal distribution for the resulting hadrons with the same slope ~ yl^e"^/^. 
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while the slope of a hadron recombining from n hard partons would steepen to 
~ B"p~^^°' (note that these n hard partons would come from n different jets!). 
On the other hand, fragmentation from a single hard parton just leads to a shift 
in the slope of the power law ~ Bp~"'^^. Given that a ~ 6 ... 8 this suggests 
that recombination of hard partons is not an important mechanism, but that 
recombination is very efficient for thermal partons. On the other hand, the 
exponential suppression of the thermal spectrum will set it in at some value of p 
and lead to a power-law spectrum of hadrons from fragmentation off jets at very 
large p, in accordance with perturbative QCD. 

Hence from very basic considerations we expect a transition from a domain 
dominated by recombination of thermal partons at intermediate px to a regime 
dominated by fragmentation of jets at very high px- It is also clear that this 
transition happens at higher values of pT for baryons compared to mesons, since 
recombination produces baryons and mesons with roughly the same abundance, 
while baryons are suppressed in jet fragmentation. 

This dual aspect of hadron production and the transition region are treated in 
different ways in the literature. 

• In publications by the FMNB group thermal recombination is supplemented 
by a perturbative calculation including jet quenching and fragmentation. 
The two components of the spectrum are simply added (|49p . No mixing 
of the thermal and hard partons is included, leading to a rather sharp 
transition between the two regions. 

• The GKL group allows coalescence between soft and hard partons as well. 
For mesons this would correspond to a term 
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• A technically very different approach is used by Hwa and Yang ()63p . Instead 

of fragmenting hard partons directly, they define the parton contents of a 
jet (initiated by a hard parton), the so-called shower distributions. They 
are given by non-perturbative splitting functions Sij{z) which describe the 
probability to find a parton of flavor j with momentum fraction z in a jet 
originating from a hard parton i. The parton content of a single jet can 
then recombine and the resulting hadron spectrum has to match the result 
from jet fragmentation. Hwa and Yang fit the shape of the parton shower 
distributions to describe the known fragmentation functions for pions, pro- 
tons and kaons (63). The power of this approach lies in the fact that the 
fragmentation part of the hadron spectrum is computed with the same for- 
malism. It is then very natural to also coalesce shower partons with thermal 
partons (j55p . 

The HY approach is also well-suited to discuss medium corrections to fragmen- 
tation in much more dilute systems like p + A collisions (64). It was found that 
the hadron-dependent part of the Cronin effect in d -|- Au collisions at RHIC can 
be attributed to coalescence of jet partons with soft partons from the underlying 
event. A more rigorous definition of parton showers and a discussion of the scale 
dependence can be found in the work by Majumder, Wang and Wang (|65p . 

2.4 Elliptic Flow 

In the momentum-space formulation it is straight forward to predict the particle- 
type dependence of the elliptic flow V2 of hadrons ()43p coming from coalescence. 
This derivation, repeated below, has been criticized as being too simplistic l\66\ 
[67|) . However, the scaling law holds numerically to very good approximation in 
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the GKL approach as weh and we will see in the next section that the data from 
RHIC follows it with surprising accuracy. We return to discussing the criticism 
further in Sec. [H 

Let us assume that the elliptic flow of a set of partons a just before hadroniza- 
tion is given by an anisotropy V2{pt) at mid-rapidity (y = 0). The phase space 
distribution of partons a can then be written in terms of the azimuthal angle (j) 
as 



where odd harmonics are vanishing due to the symmetry of the system and higher 
harmonics are neglected. / is the distribution averaged over the azimuthal angle 
(p. A general expression for the elliptic flow of hadrons coalescing from these 
partons can be derived as a function of the parton elliptic flow. For a meson with 
two valence partons a and b and for small elliptic flow t;2 <C 1 one has 



The full expressions including corrections for large elliptic flow can be found in 
()49|) . In the case of a very narrow wave function in momentum space (a oo) 
this leads to the expression 



with fixed momentum fractions Xa and Xb {xa + Xb = 1). 

Thus for hadrons consisting of light quarks which exhibit the same elliptic flow 
before hadronization we arrive at a simple scaling law with the number of valence 
quarks n: 




(13) 




(14) 




(15) 



V2{pt) = nv^{pT/n) . 



(16) 
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This scaling law had originally been suggested by several authors after first in- 
dications for scaling had been found in data gathered at RHIC ([8t[T2 |H9t [68 |) . 
Eq. (fT5]) has also been used to estimate the elliptic flow of heavy quarks from 
measurements of heavy-light systems like D mesons (I69y 60 p . The treatment has 
also been extended to harmonics beyond the second order. Generalized scaling 
laws for the 4th and 6th order harmonics have been derived in Ref. (|70p . 

2.5 Comparison of Approximations and Assumptions 

The main features and the overwhelming success of the coalescence models of 
hadronization are shared by all the approaches discussed here. However, despite 
the agreement on the general properties and their ability to describe baryon 
and meson spectra, there are different approximations and assumptions involved. 
Some have already been mentioned in the previous paragraphs. We want to 
discuss some additional points in more detail here. 

One important difference not mentioned thus far is the mass of the quarks 
in the parton phase. GKL and FMNB use effective masses that are roughly of 
the size of the constituent quark masses in the hadrons formed (i.e. rriu^d ~ 300 
MeV, nis ~ 475 MeV). This can be justified by the fact that coalescence does not 
explicitly include all the interactions. A part of the non-perturbative physics is 
encoded in the dressing of quarks, leading to a finite mass. This is also consistent 
with the requirement of (at least approximate) energy conservation. Furthermore 
quasi-particle descriptions of the thermodynamics properties of the quark gluon 
plasma estimate thermal masses of about 400 MeV (|71ll72p . However, the exact 
relation between masses in a chirally broken phase and thermal masses above 
Tc remains to be an open question. On the other hand, in the HY approach 
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massless quarks are assumed. For the phenomenology at intermediate momenta, 

p > m, masses do not play a too important role, wherefore a good description of 

measured spectra can be obtained with both assumptions. 

The missing position-space information is a weakness of the FMNB and HY 

implementations. In principle, very complex space-momentum correlations might 

exist in the parton phase before hadronization, and they might be important 

to describe elliptic flow in an appropriate fashion (|66 p 67 1[68|) . However, in the 

actual GKL computations, the spatial distribution is taken to be uniform, similar 

to the assumption used in pure momentum-space implementations. The only 

space-momentum correlation in GKL are those coming from radial flow and no 

systematic tests of more complicated space-momentum correlations are available 

in this formalism. 

On the other hand, GKL has the advantage to be able to easily accommodate 
resonance formation and decay (|53p . Direct observations of baryon anomaly and 
elliptic flow scaling are available only for stable hadrons so far. But stable hadrons 
can contain a large feed-down contribution from resonance decays, especially the 
pions ([73t [7i l[75|) . At intermediate px the role of resonance decays is somewhat 
reduced, which justifies neglecting resonances as done by FMNB and HY. The 
violation of the V2 scaling law is generally mild, which emphasizes this point. 
However, GKL shows that by including resonance decays both pT-spectra and V2 
exhibit better agreement with data towards lower pT (|76p . A schematic study 
of the elliptic flow of resonances themselves has been conducted in the FMNB 
formalism. It was found that elliptic flow is sensitive to the amount of resonances 
formed in the hadronic phase vs resonances emerging directly from hadronization 
(I77D. 
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We summarize the main differences among tiie approaclies discussed in tfiis 

section in Table [TJ 



3 Data from Elementary and Heavy-ion collisions 

In tliis section we compare various coalescence model calculations to the avail- 
able data and present predictions for future measurements. We start with single 
inclusive measurement, in particular spectra, hadron ratios and nuclear modifi- 
cation factors. We then proceed to discuss elliptic flow, particle correlations, and 
fluctuations. While all of these observables naturally focus at RHIC data taken 
during runs at ^/s = 62.4, 200 GeV, we conclude by giving an overview of the 
situation at different energies. 

3.1 Hadron Spectra and Baryon to Meson Ratios 

In Fig. Slwe show results from a coalescence model calculation of identified particle 
spectra using the FMNB method (H9]) . The spectra of neutral pions, kaons, 
protons and hyperons for central Au-|-Au collisions at 200 GeV are compared to 
date from RHIC ()38tl4ip . The salient features of the spectra are an exponential 
fall-off at intermediate pT with a transition to a harder power-law shape at higher 
Pt, as we predicted above. The transition from an exponential shape to a power- 
law shape happens at a higher pT for baryons than it does for mesons, again in 
accordance with the predictions from simple underlying principles. 

Fig.[5]shows two baryon-to-meson ratios: anti-protons vs pions (left panel), and 
A-baryons vs ^C^-mesons (right panel). Results from the GKL model for p/n^ 
(j53l) and A/2K^ (78) and from the FMNB model ^\79ij are compared to data 
from RHIC. Both calculations describe a baryon enhancement at intermediate 
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Pt that diminishes until the spectra are dominated by fragmentation at higher 

Pt- The GKL model appears to provide a better description of the data but 

a comprehensive analysis of the systematic uncertainties in the models has not 

been presented. More baryon-to-meson ratios can be found in Fig. [121 

3.2 Elliptic Flow and Quark Number Scaling 

Early identified particle measurements at RHIC showed that for px < 1 GeV/c, 
V2 at a given pj- is smaller for more massive hadrons and that when plotted vs 
niT — mo, the V2 for different species fell on a single curve. With higher statistics, 
measurements began to reveal that at higher pT the mass ordering breaks and 
more massive baryons exhibit larger V2 values ()80ll8ip . This observation led to 
the first speculation about hadron formation from coalescence and scaling of V2 
with quark number. These speculations then culminated in detailed calculations 
that we show in this subsection. 

Fig. [6] shows data on V2 scaled by the number n of valence quarks in a given 
hadron as a function of pr/n for several species of identified hadrons at ^/s^ 
= 200 GeV (|40 t l82 p . A polynomial function has been fit to the scaled values of 
V2- To investigate the quality of agreement between hadron species, the data 
from the top panel are scaled by the fitted polynomial function and plotted in 
the bottom panel. Best agreement with scaling is found for px/n > 0.6 GeV/c. 
Below that, hadron V2/n is ordered by mass. 

By combining rriT — tuq scaling and quark number scaling, one can achieve 
a better scaling across the whole momentum range (42). Fig. [7] shows V2/n vs 
(rriT — mQ)/n for several species of mesons and baryons. The scaling at low 
niT — ruQ holds with an accuracy of 5-10%. At higher niT — ruQ, a violation of the 
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simple scaling becomes apparent. In the bottom panel of Fig. [TJ V2/n is scaled 

by a polynomial fit to the meson V2/n only. The ratio of the data to the fit shows 

that baryon V2/n tends to lie below the meson V2/n. 

The break-down of the simple quark number scaling was predicted by several 

authors ([76} [66 l [67 11681 [69]) on the grounds of numerous arguments On the other 

hand, no clear consensus has emerged on whether the kinetic energy scaling 

at intermediate pT is just a consequence of pT scaling (since rwr — mo — > pt 

with increasing px) or whether it offers genuine new insights. Fig. [8] presents a 

comparison of data with the predictions for scaling violations. The bottom panel 

shows the ratio {B — M) / (B + M) , where B is V2/n for baryons and M is V2/n for 

mesons. In this figure serve as an example for mesons and A + A for baryons. 

The pion and proton V2 have been shown to be, within errors, consistent with the 

V2 of Kg and A + A, respectively ()84p83p . Two predictions for scaling violations 

are shown as well. 

Three possible sources of violations of quark number scaling have been studied 
within the GKL and FMNB implementations. One source are realistic wave func- 
tions with finite width (as opposed to the limit a = needed for the derivation 
of scaling). Both theoretical curves shown in Fig. [8] use realistic wave functions 
(note however, that in GKL in addition the quarks don't have to be collinear). 
Another correction is expected from higher Fock states which should scale with 
higher weights n + 1, n + 2, etc. A study within the FMNB framework showed 
that while thermal spectra are almost unaltered, there are visible effects for V2- 
However, those are numerically surprisingly small (jSip . Fig. [8] also contains a 
prediction including a 50% admixture of a state with one additional gluon. A 
third breaking of scaling is expected from resonance decays studied in ([76|) . 
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Those three effects cause the hadron V2/n to fah below the quark V2 values. 

The reduction is larger for baryons so that the naive scaling is broken. The 
predicted violation (|51| , l76 p are in fairly good agreement with the data ([8^. For 
> 2 — 2.5 GeV hadronization should be dominated by fragmentation, hence 
the {B — M) /{B + M) should relax to the value —0.2 if baryons and mesons from 
fragmentation have equal V2 and depend only weakly on px- The fragmentation 
contribution is not included in either theoretical calculation in Fig. [8l We discuss 
further arguments against U2-scaling in Sec. 14. 2i 

The RHIC program has also confirmed, for the first time, the existence of non- 
vanishing higher azimuthal anisotropics, beyond elliptic flow (85). The existence 
of a sizable fourth harmonic t'4 = (cos(4(/>)) had been anticipated in hydrodynamic 
calculations ((Mi)- Coalescence predictions for the relative t>4 of baryons and 
mesons provide further checks for the recombination picture. Such relations have 
been first worked out (|7Up . Concrete computations were later performed in the 
GKL model (ITS]), where it was found that the difference between baryon and 
meson is much more pronounced than for V2- This might lead to valuable 
constraints for coalescence models in future high-statistics runs at RHIC. Fits 
have been performed for identified particle from 62.4 GeV Au-|-Au collisions. 
These studies report good agreement with data for quark ^4 approximately 2x 
quark vl 

3.3 Heavy Quarks 

Coalescence has also been applied to study hadrons involving heavy quarks, in 
particular for D and B mesons ()60p 87 | l88 |) . Such studies have attracted increasing 
interest due to the surprisingly strong interaction of heavy quarks in the medium 
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first seen in the Raa (|92 U9ip and w 2 of single electrons coming from semi-leptonic 

decays of D and B mesons. While the main challenge is to understand the origin 
of this strong interaction with the medium, the hadronization mechanism plays 
a significant role in the interpretation of the data (|89 p90p . We show this in 
Figs. [9] and [10] where the Raa and the V2 of single electrons from semi-leptonic 
decays is shown together with experimental data from PHENIX and STAR l\91\ 
[92]) . Comparing the solid (coalescence plus fragmentation) and dashed band 
(fragmentation only) one notices a significant effect from coalescence. It manifests 
itself in an increase of both Raa and V2 up to pT ~ 3 GeV/c for single electrons 
(which corresponds to about pT ^ 7 GeV at the meson level. The effect is 
crucial because coalescence reverses the usual correlation between Raa and V2 
and so allows for a better agreement with the data. We note that the non- 
photonic electron spectrum can also be effected by coalescence if the Ac/D ratio 
is enhanced in Au+Au collisions compared to p+p collisions ([55]) . This is because 
the branching ratios to electrons are much smaller for charm baryons than for 
charm mesons. 

An important development is the impact of coalescence on the physics of 
quarkonia in a quark-gluon plasma. Even though coalescence has been applied to 
the J/^ for many years, the present implementations can be used to check not 
only the yield but also the spectra and the elliptic flow as a function of trans- 
verse momentum. This makes it possible to perform consistency checks between 
the spectra observed for open charm mesons and for J/^s. Such studies will be 
of particular interest at LHC where the J/^ should be dominated by regenera- 
tion in the plasma ()94[|95p . In addition, recent studies have found that even if 
the binding of a J/^ is screened in a quark-gluon plasma, the spectral function 
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still exhibits correlations above those expected for free quarks ()96p . These resid- 
ual correlations may have important implications at hadronization that can be 
studied in future recombination calculations. 

3.4 Particle Correlations and Fluctuations 

Single particle observables and elliptic flow motivated coalescence models and 
were a success story throughout the history of RHIC data taking. Later, mea- 
surements of hadron correlations challenged this picture. At RHIC it has been 
possible to measure the correlation between a trigger particle with momentum 
p^*^ and an associated particle with momentum px, typically smaller than p^*^. 
The experimental observable is usually the associated yield, which is the yield 
of correlated pairs divided by the trigger yield. Associated yields have been 
measured as a function of relative azimuthal angle Acf), and both trigger and 
associated pT (j971[98|) . This observable is ideal to detect correlations typical for 
jets. Jets give signals at Acj) = (near-side jet) and Acf) = vr (away-side jet). It 
was at first very surprising that such jet-like correlations were found with both 
trigger and associated px in the recombination domain below 4 to 6 GeV/c. 

It was then quickly realized that correlations among hadrons in this kinematic 
regime can come about through two mechanisms. First, mixed soft-hard recom- 
bination or thermal-shower recombination naturally leads to correlated hadrons. 
This was first explored by Hwa and Yang ()99p . A shower parton coalescing to 
become part of a hadron at intermediate pT will provide a correlation of this 
hadron with all those hadrons coming from fragmentation of the same jet, or the 
associated away-side jet. 

A second possibility was pointed out in a work by Fries, Miiller and Bass in 



Hadronization by Coalescence 30 
an extended FMNB framework (j621ll00tll0ip . They showed that any residual 

correlations in the tail of the bulk parton distribution automatically leads to 

correlations among the coalescing hadrons. To prove this point they introduced 

weak 2-particle correlations as corrections to the usual factorization ansatz for 

the multi-parton Wigner functions, e.g. 

Wab{Pa,Pb) = fth{Pa)fth{Pb) (1 + Cab{Pa,Pb)) (17) 

for two partons a, b. Under these specific assumptions they obtained correla- 
tions for the coalescing hadrons that are amplified by the product of valence 
quarks numbers (4, 6 and 9 for meson-meson, baryon-meson and baryon-baryon 
pairs resp.), similar to the enhancement of elliptic flow by the number of valence 
quarks. While it is not clear that the specific assumptions (very weak 2-particle 
correlations) hold at RHIC a reasonable result for associated yields as a function 
of centrality was obtained. 

Fig.[TT]shows the associated yield of near-side hadrons for trigger baryons (right 
panel) and trigger mesons (left panel) calculated in (|62|) together with PHENIX 
results from (|lU2p . A scaling law for correlations between different pairs of hadron 
species has not been observed in data so far. This is compatible with the fact 
that correlations from jet fragmentation are strong and have to be added even 
at intermediate pT, even though fragmentation is suppressed in single inclusive 
observables at the same pT (j62]) . The authors of this study argued that the phase 
space relevant for recombination at intermediate momentum is not necessarily 
completely thermalized. Rather, remnants of quenched jets, so-called hot spots 
could be an important component, leading to some residual jet-like correlations 
among partons through simple momentum conservation. Independent of the 
modeling in detail, one can conclude that recombination has been shown to be 
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compatible with measurable correlations at intermediate pT- 

Charge fluctuations (jlOSp have been shown to be consistent with the recom- 
bination process as well. They are considered to be a good probe for QGP 
formation. General expectations from coalescence are in fairly good agreement 
with data (jlOSlllO^ . A recent, more speciflc study shows that consistency with 
coalescence is obtained if the number of quarks and antiquarks is approximately 
dN/dy = 1300 for central collisions ()106p . This is in agreement with the parton 
multiplicity estimated in the GKL implementation (|53D and with the ALCOR 
model (|107p . This is a valuable consistency test for coalescence models. 

3.5 Beam Energy Dependence 

Most of the work published in the context of coalescence models focuses on 
Au+Au collisions at RHIC energies of 130 or 200 GeV. Of course, it is impor- 
tant to understand if the models can predict the correct behavior of observables, 
e.g. baryon-to-meson ratios, as a function of collision energy ^/s. Before the low 
energy Au+Au run at RHIC with = 62 GeV was completed, a prediction 
was presented within the GKL approach, utilizing a simple extrapolation of the 
model parameters (|lU8p . It was found that the p/vr ratio increases compared to 
y/s = 200 GeV while the p/vr ratio decreases. This is exactly what was measured 
when the lower energy data was analyzed (jlOOp . The predictions for scenarios 
with and without coalescence are shown together with the data in Fig. [T2l The 
data are clearly favoring the scenario with quark coalescence. The discrepancy 
found for p/vr at > 7 GeV might be due to the poor knowledge of the identifled 
proton fragmentation function (jlOSp . 

On the other side, results from heavy ion collisions at much larger energies will 
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soon become available. The Large Hadron Collider LHC at CERN will collide Pb 

ions at a center of mass energy of ^/s = 5.5 TeV per nucleon-nucleon pair. This 

will lead to QGP fireballs with much higher temperatures. On the other hand one 

can also predict that the number of hard processes increases tremendously due to 

the rising gluon distribution at small Bjorken-x. Naively, one would expect the 

window in pT where coalescence is dominating to increase. However, the estimates 

for this region depend delicately on the radial flow (which pushes the coalescing 

hadrons to higher px) and jet quenching (which leads to less fragmented hadrons 

at high pt)- 

Possible scenarios have been explored in the FMNB framework using different 
assumptions for the radial flow (|110p . These estimates are shown in Fig. [T3j 
Recently, a more systematic study of elliptic flow as a function of collision energy 
was published 

4 Challenges and Outlook 

The RHIC program has provided remarkable evidence that coalescence of quarks 
is the dominating mechanism for hadronization from a deconfined plasma. Nonethe- 
less, some problems remain unsolved and several new questions are raised by the 
formalism itself. E.g. it is appealing to apply recombination at low momenta 
where the phase space is more dense. Some of these problems have been touched 
upon briefly in previous sections. We will discuss them in more detail below. 

4.1 Energy and Entropy 

A basic issue that involves all approaches based on instantaneous projection is 
that of energy conservation. The underlying kinematics of the projection is ef- 
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fectively 2 — > 1 and 3 — > 1, which makes it impossible to conserve 4- momentum. 

This is somewhat mediated at intermediate transverse momenta, pT > fn, where 
the kinematics is essentially collinear and violations of energy conservation are 
suppressed by factors m/pT or kT/pT where kT is the intrinsic transverse mo- 
mentum of a parton inside the hadron. In principle, this is not really acceptable, 
since the formalism should be easily extendable to low px where coUinearity is 
missing. In fact, a smooth matching with bulk coalescence models like ALCOR 
should be possible, which describe multiplicities and related observables at low 
Pt successfully ()30U112p . Interestingly, a naive extension of the GKL approach to 
low momenta does not lead to striking disagreement with the experimental data 
(|53 U 89 p . However from the theoretical point of view the issue of imperfect energy 
conservation is clearly unsatisfying. 

Energy conservation has to be achieved through interactions with the surround- 
ing medium. Naturally, approximations to this multi-particle dynamics have to 
be applied to make the problem tractable. One way is to introduce an effective 
mass distribution for the quarks as a way to incorporate some in-medium effects 
(jll3p . This allows to enforce both momentum and energy conservation and one 
finds fairly good agreement with data for p^^-spectra. 

A promising new and very powerful approach has recently been developed by 
Ravagli and Rapp (RR) ([38|) . They replace the instantaneous projection of quark 
states onto hadron states by a procedure which solves the Boltzmann equation 
for an ensemble of quarks which are allowed to scatter through hadronic states. 
Thus the hadrons are given through cross sections with a certain width. This 
implementation naturally conserves 4-momentum. Ravagli and Rapp find quite 
good agreement with data for p^^-spectra. They also confirm V2 scaling (neglecting 
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position-momentum correlations as the other approaches). However, they find 

that kinetic energy scaUng {v2 vs rriT — Jtt-o) is in even better agreement with 

experimental data, cf. Fig. [71 The RR formalism with energy conservation is the 

only one really suited to address the question of kinetic energy scaling. 

A related issues is entropy conservation. Coalescence through instantaneous 
projection seems to reduce the number of particles by about a factor two, which 
understandably rises the question whether the second law of thermodynamics is 
violated. However, strictly speaking this formalism should only be applied at 
intermediate px where only a small fraction of the total particle number (< 2%) 
is located. Furthermore, the situation is much less dire if resonance production 
is taken into account which significantly increases the number of hadrons in the 
final state ([S ^ISS]) . Entropy depends not only on the number of particles, but 
also on the degeneracies in both phases and on the masses. 

In addition one should also take into account the interaction among quarks. It 
has been shown for an isentropically expanding fireball, using the lattice equation 
of state, that the evolution of the effective number of particles reduces signifi- 
cantly around the crossover temperature (jll4p . This could help to solve the 
entropy problem inherent to instantaneous quark coalescence, as also pointed 
out by Nonaka and collaborators ()106p . However, it is still a challenge to find 
a consistent approach to conserve both energy and conserve or increase entropy, 
together with a good description of single particle spectra and elliptic fiow for 
both low and intermediate pT- 
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4.2 Space-Momentum Correlations 

An important open question is the relation between space-momentum correlations 
and V2 scaling. The valence quark number scaling of elliptic flow was derived in 
a pure momentum-space picture. This means that the scaling has been explicitly 
proven only if the coalescence probability is homogeneous in space. GKL have 
gone one step beyond by including correlations of radial flow with the radial 
coordinate r. They find that scaling still holds to a good approximation with 
some small violations (j53y 76 1[89|) . 

However, the situation could be very different if more realistic correlations of 
flow with the spatial azimuthal angle (p are taken into account. One should expect 
a strong correlation between the spatial azimuthal angle if and the momentum 
azimuth (f). A detailed discussion of effects coming from space-momentum corre- 
lations can be found in the work by Pratt and Pal (66). They also map out a 
class of phase space distributions that lead to approximate scaling. 

Parton cascade studies that calculate the time evolution of the phase space 
distributions find that approximate scaling between baryons and mesons still 
persists even if strong deviations of V2 at the quark level are seen ([57| . However 
it is not clear how this depends on the freeze-out criteria, on the width of the 
wave functions and on the interplay with jet fragmentation. Another study on 
the effect of phase space distribution can be found in Ref. (|68p . 

Small violations of V2 scaling have been observed, but as discussed in detail in 
Sec. [3] they can be explained solely by wave function effects, resonance contribu- 
tions and contributions from higher Fock states in hadrons ()5Hll081lll5p . If the 
scaling feature were accidental and strongly dependent on details of the phase 
space distribution, the very different dynamical evolution at LHC might lead to 
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much stronger scaling violations there. A better understanding of HBT measure- 
ments could also supply fundamental information on this issue. It would also be 
interesting to see how realistic phase space correlations will fare in a dynamical 
coalescence model like the RR formalism. 

4.3 Outlook 

Quark coalescence models for heavy ion collisions have reached a certain level of 
maturity, but it has also become clear that there are limitations. We hope that 
several issues will attract attention in the future. 

Within the established projection formalism several open questions can be 
addressed. A huge amount of data on 2- and 3-hadron correlations has been 
collected. While preliminary studies have shown that correlations are in prin- 
ciple compatible with recombination, a comprehensive effort to understand the 
data in a picture that contains jets, jet quenching, and coalescing partons at 
intermediate pT still has to be developed. It would have to include a realistic 
microscopic modeling of the coupling between the medium and jets and how jet- 
like correlations can be conferred to the medium. A second issue concerns the 
role of resonance production. Little is know about the relative probabilities of 
coalescence into stable hadrons and unstable resonances. As we have seen above 
this is an important issue for multiplicities and entropy production as well as V2 
scaling violations (in particular for pions). 

Dynamical transport implementations like the one developed by Ravagli and 
Rapp are very promising candidates to investigate more fundamental open ques- 
tions. E.g., it would be straight forward to implement resonances and stable 
hadrons in a realistic fashion. Progress could be made on the issues of kinetic 
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energy vs transverse momentum scaling of V2, the role of space-momentum cor- 
relations for elliptic flow scaling, and entropy production. There is also a need 
to explore dynamical coalescence coupled to realistic transport models for the 
parton and hadron phase. 

There is a list of more profound questions which we have not touched upon 
yet at all. Coalescence of particles can be found in systems which do not exhibit 
confinement (e.g. in plasmas of electrons and protons). Confinement does not 
play a big role in any of the current implementations of coalescence. (In parton 
cascades, non-coalescing partons are usually fragmented, the only tribute to the 
fact that there are no free partons allowed in the vacuum.) Nevertheless there 
should be a fundamental difference between confining and non-confining theories. 
Transport implementations need to explore this difference in the future. 

It is also not clear what the role of chiral symmetry breaking during the co- 
alescence process is. Most implementations give constituent-like masses to the 
quarks, but no direct connection to chiral or thermal masses is made. Unfortu- 
nately the current observables do not seem to be sensitive to the nature of the 
quark masses. We would hope that improved implementations together with new 
high-statistics data might allow us to address this question. 

5 Conclusions 

The first stage of the RHIC program has provided clear evidence that hadroniza- 
tion at transverse momenta of several GeV/c is modified when compared to p+p 
collisions in the light quark sector. The available data is only compatible with a 
hadronization process through coalescence of quarks. The baryon enhancement 
and the robust scaling of the elliptic flow with the number of valence quarks are 
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signatures which rule out other explanations. 

We have presented a comprehensive overview of the available coalescence mod- 
els, which are mostly based on an instantaneous projection of quark states onto 
hadron states. On the other hand, dynamical coalescence uses scattering of quark 
into hadron states in a transport approach. We have discussed some of the weak- 
nesses of current implementations and how the field might evolve in the future. 

Let us conclude by discussing a statement that has naturally arisen after coa- 
lescence models had been applied to RHIC. It was pointed out again and again 
that coalescence might be the most convincing argument to show that confine- 
ment takes place at RHIC and a quark gluon plasma is indeed formed. The 
argument used relies on the fact that elliptic flow V2 is a collective effect (coming 
from the hydrodynamic expansion due to pressure gradients), and that this col- 
lectivity seems to happen on the parton level, leading to a universal elliptic flow 
for quarks just above Tc. In other words, elliptic flow of hadrons at intermediate 
Pt did not emerge from hadronic interactions. 

This is indeed very remarkable and it is a strong argument for deconfinement. 
All signatures for deconfinement use indirect arguments and need some kind of 
theoretical input to reach this conclusion. Coalescence, and in particular the V2 
scaling, appear to be convincing because almost no additional assumptions seem 
to be needed. We hope that this argument is solidified with future improvements 
in our understanding of data and of the mechanism of recombination. 
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Figure 1: Left: p/vr^ ratios measured in central Au+Au collisions at ^/s^ = 
200 GeV at RHIC, compared to measurements from + e~ and p + p collisions. 
Right: The ratio A/2Kg for central and mid-central Au+Au collisions at ^/s^ = 
200 GeV measured by STAR. The p/tt~ ratio from p+p collisions from STAR is 
shown for comparison. 
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Table 1: A summary of key differences between the most popular implementa- 
tions: GKL = Greco, Ko, Levai; FMNB = Fries, Miiller, Nonaka, Bass; HY = 
Hwa, Yang. 
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Figure 2: Nuclear modification factors {Rcp) for various identified particles mea- 
sured in Au+Au collisions at ^/s^ = 200 GeV by the STAR collaboration. The 
and A + A Rcp values demonstrate that strange baryon yields are enhanced 
in central Au+Au collisions compared to strange meson yields. Later, measure- 
ments of the i;^, S + S and i7 + 17 showed that the rate of increase of the particle 
yields with collision centrality depended strongly on whether the particle was a 
baryon or meson with the mass dependence being sub-dominant: the baryon and 
meson Rcp values fall into two separate bands (indicated by lines to guide the 
eye) with the baryon Rcp larger than the meson Rcp- 
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Figure 3: V2 for a variety of particles from a minimum-bias sample of Au+Au 
collisions at ^/s^ = 200 GeV measured by the STAR (6) and PHENIX ([35]) col- 
laborations. Curves show the results from hydrodynamic model calculations (jl4p . 
f 2 values also show that baryon production at intermediate pT is enhanced in the 
in-plane direction, leading to larger baryon V2- This observation is incompatible 
with expectation of V2 coming from parton energy loss. 
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Figure 4: Hadron pT-spectra at midrapidity from 200 GeV central Au+Au 
collisions. The curves show the recombination and fragmentation components of 
the spectra obtained in the FMNB formalism along with the total which compares 
well with the data. 
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Figure 5: Ratios of baryon yields to meson yields for central Au+Au collisions 
at 200 GeV. The GKL and FMNB calculations for p/tt~ (left) and A/2K^ are 
compared to STAR and PHENIX data. 
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Figure 6: Top panel: The elliptic anisotropy parameter V2 scaled by quark 
number n and plotted vs pt/^. A polynomial curve is fit to all the data. The 
ratio ol V2/n to the fit function is shown in the bottom panel. 
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Figure 7: Quark number scaled elliptic flow vs {rriT — "m-o) /n. In the low — tjiq 
region, the scaling is improved by plotting vs tjit — ruQ. All data is fit by a 
polynomial curve and the ratio V2/n to the fit function is shown in the bottom 
panel. 
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Figure 8: Top panel: Quark number scaled V2 showing violation of ideal scal- 
ing. A polynomial is fit to all the available data. Bottom panel: The difference 
between quark number scaled baryon V2 and quark number scaled meson V2 di- 
vided by the sum: {B — M)/{B + M). The ratio is formed using hyperons and 
kaons. The solid curve shows model predictions from FMNB using realistic wave 
functions and a 50% admixture of a higher Fock state containing an additional 
gluon. The dashed line shows calculations in the GKL model with realistic wave 
functions. 
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Figure 9: Nuclear modification factor Raa of single electrons from semi-leptonic 
decays in Au+Au collisions at 200 GeV. The solid line represents the predictions 
from a coalescence plus fragmentation model ()89p for electrons from D and B 
mesons (shaded bands) and from D mesons only (lines). The shaded band reflects 
the theoretical uncertainty in the heavy quark diffusion coefficients (87). The 
dashed lines are the results without coalescence. The data are taken from (|9ip92p . 
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Figure 10: Elliptic flow V2 of single electrons from semi-leptonic decays in Au+Au 
collisions at 200 GeV. The lines represent the same calculations as in Fig. [9l The 
data are taken from Ref. (j9ip 
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Figure 11: Associated hadron yields on the near-side as a function of number of 
participants for meson triggers (left) and baryon triggers (right) from (j62p . The 
diamonds represent the expected hadron correlations if fragmentation is the only 
source of correlations and recombination is correlation-free. Triangles show the 



same calculation with small 2-particle correlations among coalescing partons. 
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Figure 12: p/vr^ ratio and p/tt" ratio in central Au+Au collisions at 62.4 GeV. 
The predictions of the GKL model (coalescence plus fragmentation) (jlOSp are 
shown by thick solid lines for p/ir^ and by thick dashed lines for p/ir^; the 
prediction from fragmentation only are the corresponding thinner lines. The 
data from STAR are taken from Ref. ()109p 
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Figure 13: Predictions for vr*^ (left) and proton (right) spectra in central Pb+Pb 
collisions at LHC using radial flow parameters (3 = 0.65,0.75 and 0.85 resp. The 
larger the radial flow the more the recombination region extends to higher px, 
possibly up to 10 GeV/c (from (fTTOl) ). 



